Global Subsonic and Subsonic-Sonic Flows 
through Infinitely Long Axially Symmetric Nozzles 



Chunjing Xie* and Zhouping Xin^ 

Abstract: In this paper, we establish existence of global subsonic and subsonic-sonic 
flows through infinitely long axially symmetric nozzles by combining variational method, 
various elliptic estimates and a compensated compactness method. More precisely, it is 
shown that there exist global subsonic flows in nozzles for incoming mass flux less than 
a critical value; moreover, uniformly subsonic flows always approach to uniform flows at 
far fields when nozzle boundaries tend to be flat at far fields, and flow angles for axially 
symmetric flows are uniformly bounded away from ir/2; finally, when the incoming mass 
flux tends to the critical value, subsonic-sonic flows exist globally in nozzles in the weak 
sense by using angle estimate in conjunction with a compensated compactness framework. 

Keywords: Axially symmetric nozzles, Subsonic flows, Subsonic-sonic flows, Stream 
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1 Introduction and Main Results 

This is a continuation to [13J on the study for subsonic and subsonic-sonic flows in multi- 
dimensional nozzles. Two dimensional subsonic and subsonic-sonic flows through infinitely 
long nozzles were studied in detail in [13J. In particular, global smooth subsonic flows and 
their properties have been established for incoming mass fluxes less than the critical value, 
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while existence of subsonic-sonic flows is proved for the critical mass flux. However, many 
of arguments are special to 2-dimensional flows, and it seems difficult to generalize them to 
the more realistic 3-dimensional flows. 

In this paper, we would like to investigate the 3-D flows in nozzles which are infinitely 
long and axially symmetric. 

As far as axially symmetric flows are concerned, one should note the significant result 
due to Gilbarg, [7\, where he showed that if an axially symmetric subsonic nozzle flow 
approximates to uniform flows at far fields, then the flow speed on the boundary is monotone 
increasing with respect to the incoming mass flux by a comparison principle, however, 
existence of such flows is not known. For free boundary problems, in [I], Alt, Caffarelli and 
Friedman gave a complete study for flows with jet and cavitation by variational methods. 

Let us start with 3-D isentropic compressible Euler equations, 

(Pu)x + (pv)y + (pw) z = 0, 

(pu 2 ) x + (puv)y + (puw) z + p{p) x = 0, 

(1) 

{puv) x + (pv 2 ) y + (pvw) z +p(p) y = 0, 

{puw) x + (pvw)y + (pw 2 ) z + p{p) z = 0, 
where p is the density, (u, v, w) is the velocity, and p = p(p) denotes the pressure. In 
general, it is assumed that p'(p) > for p > and p"(p) > 0, where c(p) = \/p'(p) is called 
the sound speed. Important examples include polytropic gases and isothermal gases, for 
polytropic gases, p = Ap 1 where A is a constant and 7 is the adiabatic constant with 7 > 1; 
and for isothermal gases, p = c 2 p with constant sound speed c. 
Suppose that the flow is also irrotational, i.e. [5], 

Uy -v x = 0, u z - w x , v z -w y = 0. (2) 

Then it follows from ([1]) and ([2]) that the flow satisfies Bernoulli's law 

where q = Vu 2 + v 2 + w 2 , and C is a constant depending on the flow. There are some 
basic facts about irrotational isentropic steady flows, see [5], which are consequences of 
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Bernoulli's law (|3|). First, p is a decreasing function of q, attains its maximum at q = 0. 
Second, there is a critical speed q c such that q < c(subsonic) if and only if q < q c . Finally, 
pq is a nonnegative function of q, for q > 0, which is increasing for q £ (0, q c ) and decreasing 
for q > q c , and vanishes at q = 0. so attains its maximum at q = q c , therefore, that the 
flow is subsonic is equivalent to pq < p c q c and p > p c . Therefore, we can nondimensionalize 
the flow as in jSJ [13], such that q cr = 1, p cr = 1, and then Bernoulli's law ([3]) reduces to 

2 Ji p 2 
Since p'(p)/p > for p > 0, so (JH) yields a representation of the density 

p = g(q 2 ), (5) 

moreover, g is a decreasing function. For example, for polytropic gases, after the nondi- 
mensionalization, p = /9 7 /7, and ([5]) is nothing but 

2 ( 1 + l-( 1 -l)q 2 \^ 

p = g(q) = [ £ I • (6) 

Furthermore, p is a two- valued function of (pq) 2 - Subsonic flows correspond to the 
branch where p > 1 if (pq) 2 G [0, 1). Set 

P = H((pq) 2 ) (7) 

such that p > 1 if (/?(/) 2 £ [0,1), therefore, -ff is a positive decreasing function defined on 
[0,1], twice differentiable on [0,1), and satisfies H(l) = 1. Moreover, it follows from ([5]) 
and ([7]) that (pq) 2 is given in terms of q 2 as 

(pq) 2 = G(q 2 ). (8) 

Thus 

g(q 2 ) = H(G(q 2 )). 
Suppose that the wall of nozzle is impermeable so that 

(u, v, uu) ■ n = 0, (9) 
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where ~n is the outward normal of the solid boundary. 

Due to ([2]), one can introduce a velocity potential for the flow such that 

$ x =u, $ y = v, <f> z = w. 

Thus the continuity equation becomes 

div( 5 (|V$| 2 )V$) = 0. 

Assume now that the nozzle is axi-symmetric as given by 



D = {(x,y,z)\0 < \/y 2 + z 2 < f(x), -oo < x < oo}. (10) 

Consider a smooth flow in the nozzle. Then it follows from continuity equation and (|9|) that 
mass fluxes through each section which is transversal to the symmetry axis are the same. 

Thus, the problem of finding solutions to smooth flows in a 3-D nozzle reduces to solving 
the following problem, 

div( 5 (|V$| 2 )V$) = 0, inD, 

§1 = 0, on dD, (11) 

k f s9 (\V^)m,dS = m 0l 

where S is the surface transversal to the axis, and I is the normal to S which directs to the 
positive axial direction. 

In this paper, it is assumed that there exists a £ (0, 1) such that 

H/'llc7«(M) < oo, and inf / = b > 0. (12) 

Now the main results of this paper can be stated as follows. 

Theorem 1 Suppose that the nozzle boundary satisfies Uty). Then there exists a positive 
constant fh depending only on f such that if mo < m, there exists an axially symmetric 
uniformly subsonic flow through the nozzle. More precisely, there exists a smooth solution 
$ G C°°(D) to (TUP such that 

sup|V$|<l, (13) 

D 



and 



u(x,y,z) = ® x = U(x,r), v(x,y,z) = $ y = V(x,r)-, w(x,y,z) = $ z = V(x,r)-, (14) 



where r = \/y 2 + z 2 and U(x,r), V(x,r) are smooth in their arguments, and V(x,r) van- 
ishes on the symmetry axis. 

If the wall of the nozzle tends to be flat at far fields, for example, rescaling if necessary, 
one may assume that 



f(x) — ► 1 as x — > — oo, fix) -^a>0asi;-> oo. 



(15) 



Then the following sharper results hold. 



Theorem 2 Suppose that the wall of the nozzle satisfies both lllty) and M5\l. Then there 
exists m > such that if < uiq < m, there exists a unique axially symmetric uniformly 
subsonic flow through the nozzle with the properties that 



M(m ) = sup |V<J>| < 1, 



and 



(U,V)-({G-H^)} 1/2 ,0) 



7T Z 



(U,V) - ({G-\^L)}^\0) 



as x — > — oo, 
-> as x — ► oo, 



(16) 



uniformly in r, where G is defined by |2]); moreover, M(mo) ranges over [0, 1) as uiq varies 
in [0, m). Furthermore, ifO< rriQ < m, i/ie axia/ velocity is always positive in D, i.e. 



u > 0, 



(17) 



and, the flow angle, uj = arctan jj , satisfies 



U) < to < U), 



(18) 



where 



u_ = minjinf arctan f'(x), 0}, to = maxjsup arctan f'(x), 0}. 



(19) 



Moreover, for any given m £ (0,m), there exist a positive constant 5 = 5(m) > 0, such that 
if m £ [m, rh), then 

q(m) = inf |V*| > 5. (20) 

We now study the limiting behavior of these subsonic flows in the nozzle when then 
the cross-section mass fluxes mo approaches the critical value. In fact, as mo f m, the 
corresponding flow fields tend a limit which yields a subsonic-sonic flow in the nozzle. 

Theorem 3 Let {mo, n } be any sequence such that mon — > rh as n — > +00. Denote by 
(U n , V n ) the global uniformly subsonic flow corresponding to mo, n o,s guaranteed by Theorem 
2. Then there exists a subsequence, still labelled by {(U n ,V n )} associated with {mo, n } such 
that 

U n ^U, V n -> V, (21) 
9{ql)Un -» 9(q 2 )U, g{q 2 n )V n -» g(q 2 )V, (22) 

where q n = U 2 + V%, q 2 = U 2 + V 2 , and g(q 2 ) is the function defined by f5|] through 
Bernoulli's law, all the above convergence are almost convergence. Moreover, this limit 
yields a three dimensional flow with density p(x,y,z) = g(q 2 )(x,r) and velocity 

y z 
u(x,y,z) = U(x,r), v(x,y,z) = V(x,r)-, w(x,y,z) = V(x,r)-, 

where r = \/y 2 + z 2 , which satisfies 

u y — v x = 0, v z — w y = 0, w x — u z = in D 

in the sense of distribution, moreover, for any n £ C^°(D) 

j j j (pu, pv, pw) ■ Vrjdxdydz = 0. 

This implies that (u, v, w) satisfies boundary condition (0J) as the normal trace of the diver- 
gence field (pu, pv, pw) on the boundary. 

Before we prove the theorems, there are a few remarks in order. 
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Remark 1 In contrast to two dimensional plane flows, three dimensional flows are much 
more complicated. Indeed, some of the key arguments in [13] can not be applied to three 
dimensional case directly. Even for irrotational steady axially symmetric subsonic flows, 
there are some difficulties near the symmetry axis, see (|29p . Therefore, it seems difficult 
to show the existence of subsonic flows by fixed point argument as in plane flows in |13j . 
Fortunately, for axisymmetric flows, equation (|29p has a variational structure, which is one 
of the key points to show the existence of subsonic solutions. 

Remark 2 It should be noted that one cannot adapt the analysis of [1] directly to study 
the properties of the subsonic flow in Theorem [2] since for jet flow, the pressure is prescribed 
on the jet surface, so the flow speed is known by Bernoulli's law, thus it is easier to see 
whether the flow is subsonic and whether it approaches to uniform flows at far fields. 

Remark 3 In all the theorems in this chapter, we require only C 1,a smoothness of /. 
Similar to the proofs given in this paper, one can prove all results in [13] under the condition 
that nozzle boundaries are C 1,a instead of C?'°. Furthermore, it is only required that / 
itself tends to constants at far fields instead of its higher derivatives, which improves the 
results in [T3] . 

Remark 4 Theorem [2] provides the existence of flows studied by Gilbarg in [7]. Moreover, 
applying the comparison principle obtained by Gilbarg in [7], the maximum speed of flows 
obtained in Theorem [2] is monotone increasing with respect to incoming mass flux. 

Remark 5 There are some fragmentary descriptions of some phenomena on the axially 
symmetric subsonic flows past a body, for the reference, please refer to [3], [9], [8]. For 
applications of the theory of compensated compactness to two dimensional transonic and 
subsonic-sonic flows, please see [12], [4], [13] . 

The rest of the paper is organized as follows: in Section 2, we derive the governing 
equation and boundary conditions for axially symmetric irrotational flows. In Section 3, 
we adapt the variational method used in pQ to prove Theorem 1. Subsequently, in Section 
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4, we prove that subsonic flows will approach uniform flows at far fields when the nozzle 
boundaries tend to be flat at far fields, which will yield the existence of the critical value 
for incoming mass fluxes. In Section 5, positivity of axial velocity and uniform estimates 
for flow angles for axially symmetric flows are proved. In last section, Section 6, we use a 
compensated compactness framework to show the existence of weak subsonic-sonic flows. 

2 Axially Symmetric Flows 

In this section, we will derive the governing equations and boundary conditions for axially 
symmetric irrotational flows in cylindrical coordinates and in terms of stream functions. 

In the cylindrical coordinates (x, r, 9), let the fluid density and velocity be p(x, r, 9) and 
(U(x, r, 9), V(x, r, 9), W(x, r, 9)), where U, V, and W are axial velocity, radial velocity and 
swirl velocity respectively. Then (x,y,z), p, and (u,v,w) satisfy 

x = x, y = r cos 9, z = r sin 9; 
P(x, y, z) = p(x, r, 9), u(x, y, z) = U{x, r, 9); 
v(x, y, z) = V(x, r, 9) cos 9 + W(x, r, 9)(— sin 9), 
w(x, y, z) = V(x, r, 9) sin 9 + W(x, r, 9) cos 9. 

It should be noted that for axi-symmetric flows, one has 

U(x,r,9) = U{x,r), V(x, r, 9) = V(x, r), W(x, r, 9) = W{x, r). 

Since the flow is also assumed to be irrotational, one has 

(rW) r 

V Z - Wy = = 0, 

r 

this implies that 




Thus W = since W is bounded near r = 0. Therefore, for axially symmetric irrotational 
flows, one has 

u = U(x,r), v = V(x,r)-, w = V(x,r)-, and p = p(x, r), (23) 
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where r = \Jy 2 + z 2 . Then the continuity equation reduces to 



(rpU) x + (rpV) r = 0. 



Moreover, the irrotational condition (|2|) changes to 



U r - V, = 0. 



Bernoulli's law is still of the same form with q = y/U 2 + V 2 . 

Due to (I24p . one can introduce a stream function ip = tp(x,r) such that 

ipr = rpU, 4> x = -rpV. 



Then Bernoulli's law ([4j) becomes to 
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'i P 2 
Therefore, it follows from (JTj) that p can be represented as 



p = H 



so the irrotationality (|25p changes to 



div \H 



The no-flow boundary condition ([9]) on the nozzle wall becomes 



(U, V) ■ N = 0, 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 



(30) 



where N is the normal of the curve r = f(x). It follows from (|30p that ^ is a constant in 
each connected component of the solid boundaries. 

Note that for smooth axisymmetric flows in the nozzle, it follows from (]26p that ip is a 
constant on the symmetry axis. Thus r = is a streamline. 

Since the flow is axially symmetric, one may consider only symmetric part of the domain. 

Let 

Q, = {(x,r)\0 < r < f(x), -oo < x < oo} (31) 



with boundaries 



T\ = {(x,r)|r = 0, — oo < x < 00} , T 2 = {(x,r)\r = f(x), — 00 < x < 00}. (32) 



For convenience, we denote by Dq the three dimensional domain induced by f2, 



Dq = {(x, y, z)\0 < a/ 'y 2 + z 2 < f(x), —00 < x < 00}. 



(33) 



Then, to study the 3-dimensional problem, (jlip . for axisymmetric flows, one may first 
study the following 2-dimensional problem 



div \ [H 



r 



^ = 0, 



m = 0, in n, 

on Ti, 
on T 2 . 



(34) 



3 Subsonic Flows Associated with Small Incoming Mass Flux 

This section is mainly devoted to the proof of Theorem 1. Our approach is motivated 
strongly by the important work pQ by Alt, Caffarelli and Friedman. The proof can be 
divided into 10 steps. 

Step 1. Subsonic truncation and shielding singularity. By direct calculations, it is easy 
to find that the derivative of function H(s) goes to negative infinity as s — > 1. To control 
the ellipticity and avoid singularity of H' , one may truncate H as follows 



His) 



H(s), 



if < s < m 2 , 



(35) 



where rh < 1 is a constant to be determined, and H is a smooth decreasing function. Set 



q 2 = s/H 2 {s). 



(36) 



Since H 2 {s) — 2HH'{s)s > 0, we can represent s as a function of q 2 , s = G{q 2 ). Obviously, 
G is an increasing function. Define p = g(q 2 ) as 



~g{q 2 ) = H{G(q 2 )). 
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(37) 



Then it is easy to check that 



H(G(q 2 )) 



> v 



(38) 



dq 2 H(G(q 2 )) - 2H'(G(q 2 ))G(q 2 ) 

for some positive real numbers A and v > which depend on H. 

To treat the singularity in the coefficients of the equation (|29p as r — ► 0, one may shield 
the singularity by first solving the following problem 



div H 



r+S 



r+S 



0. 



in f2. 



(39) 



ijj = 0, on T%, 

ip = m, on T%. 

Step 2. Variational problem. The problem (|39p is a boundary value problem for an 
elliptic equation in a unbounded domain, therefore, we use a series of Dirichlet problems in 
bounded domains to approximate it. Thus consider first the following problem 



div 

tf,= 



H 



Vi/> 



r+5 



r+6 



o. 



m. 



in Ql, 
on drt^, 



(40) 



where £Il = {{x, r)\(x, r) E SI, \x\ < L}. The problem, (|40p . can be solved by a variational 
method. The existence of solution to problem (|40p is equivalent to find minimizer ijj 5 L E 



= {</)|(/) E W 1,2 (fiz,), — itjt^ £ ' 2 (^l)} f° r the following minimization problem 



inf M<t>), 
<t>&A L 



where 



and F is defined by 



F 



r + J 



(r + 5)dxdr, 



F(s) 



Since H is a smooth decreasing function, therefore, J r (pi,P2, r ) = F 



(P1.P2) 

r+<5 



(41) 



(42) 



(43) 



(r + 5) is 



a convex function of p = (pi,P2), using the standard theory in calculus of variations, for 
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example, Theorem 2 and Theorem 3 in Section 8.2 in [6], the problem (|4ip has a unique 
solution since the functional Jl is also coercive. 

Step 3. Estimates for minimizers. For each L, there exists a unique solution ij) L of 
the problem (|4ip . Each minimizer ipi is a weak solution to the Dirichlet problem of the 
Euler-Lagrange equation, (|40p . Then by a weak maximum principle for the problem (|40p . 
see Theorem 8.1 in [10] . one gets 

< ip S L < rn in Vt L . (44) 

Using Caccioppoli's inequality, both in interior and on the boundary, and Theorem 6.5 and 
Theorem 6.8 in [11], one can obtain 

l|V^||r» ( n Ji )<C(L ) [|/ , [| c o,<5 > m) I Vfe > 2L. (45) 



Then Holer estimates for the gradient of minimizers to the functional (|42p . and Theorem 
8.6 in [TT] . imply that there exists a\ G (0,a) such that 

Ui\\c^{n L )<C{L : a u \\f\\c^5,m), Vk > 2L. (46) 

Moreover, the interior Schauder estimate, Theorem 10.18 in [TT], shows that for any X CC 
Ol, it holds that 

\\1>k\\c*'"ip) <C(E,L, ai ,S,m), Vk>L. (47) 

To recover the singularity later by taking the limit 5 — > 0+, we need a more precise 
estimate than (|4*4"j) . Set 

- (r + 5) 2 

= -^P-™' ( 48 ) 

where b is defined in (|12|) . Then it is easy to check that ^ satisfies the equation 

-l 



div \H 



r + <5 



r + 6 



0. (49) 



Because fi^ satisfies a uniform exterior cone condition, tfi^ G C°(r2i) by Theorem 8.29 in 
[TO] . Moreover, by (|4T|) . V'i S C 2,ai (f2L). Therefore, both and -0 satisfy the equation 
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in (|40p on in the classical sense. Obviously, ip > ip s L on (90^. Thus, it follows from a 
comparison principle, Theorem 10.1 in [10], that 



rp S L <ip inn L . (50) 

Step 4. Existence of solutions to (|39p . By a diagonal process and Arzela-Ascoli lemma, 
it follows from (|46[) that there exists a sequence {n^} such that 

<Xn„ fc - / m C 1 '^) for VL > 

with < \i < ai. Therefore, ^ 5 is a weak solution to the problem (|39p . Then it follows 
from g7]) that ij) S G C 2 ^(ft L ), VL > 0. 

Step 5. Recover singular coefficients. Due to (f5T)j) . we have 

< S( r + 5 ) 2 - 

cr 

Therefore, ^ 5 ^ on Ti. Moreover, for Ve > 0, on each set Ql,e = { {x, r)\ \x\ < L,e < r < f(a 
it follows from Caccioppoli's inequality and Holder gradient estimate in a similar way as for 
15"]) and (|46p . that V' 5 satisfies the following estimate 



U S \\o^i(0^)< C ( L ^A\f'\\c^rri). (51) 

Due to a diagonal process and Arzela-Ascoli Lemma again, there exists a subsequence {5k} 
such that 

iP 5k -> ^ in C 1,M (Oi, £ ) for each L > 0, e > 0. 

In particular, 

■i/A — > ?/> pointwise in S7. (52) 
Moreover, tp £ C 1 '^ solves the problem 



div [H 



^ = 0, in a 



(53) 



2 



weakly and satisfies 



< V < jir 2 . (54) 
o 2 
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It follows from the standard bootstrap arguments that tp G C 2,tl (Q). 

Step 6. Subsonic estimate near the symmetry axis. In this step, our aim is to show that 



(x, r) 

r 



< Cm 



(55) 



for < r < §. To do this, we note an important observation due to [T], that if r$ < 6/2, 
then 



ipo{x,r) = pil)(x +tx,r Q + tr), t- } . 



(56) 



satisfies 



div ^ l^H 

It follows from (|54|) that 



V^o 



2 + r 



2 + r 



in Si ((0,0)) 



0<Vo< Cm in £i ((0,0)). 



Therefore, by Moser's iteration, Theorem 8.18 in [10J, one can get 



< Cm in Si ((0,0)). 



In particular, 

^(x ,r ) = |VVo(0,0)| < Cm. (57) 
r 

Step 7. Subsonic estimate away from the symmetry axis. In this step, we derive the 
estimate (f55|) for r > 6/4. For any given (£cb r o) G ^oo,b/4 = {(^j r )l( x > r ) G > 6/4}, 
noting that all the coefficients in the equation ([53]) are bounded on B b / s ((xo, ro)) C\ fio^&M, 
moreover, 

< ip < m 

due to (|4"4"|) and ([52]) . one can derive by Moser's iteration that 

|Vy>(x ,ro)| < Cm. 



Therefore, 



since ro > 6/4. 



(zo,?U 



< Cm, 



(58) 
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Step 8. Uniform Holder continuity of velocity field near the symmetry axis. It follows 
from Step 6 and Step 7 that the flow is subsonic except on the axis when the incoming 
mass flux is sufficiently small. To show that the flow is subsonic globally, we first need to 
show that the velocity field is well-defined along the symmetry axis. In fact, we have the 
following stronger results. 



Lemma 4 Let if) be a solution to problem 153\) satisfying 



<C inn. 



Then — is uniformly Holder continuous up to the symmetry axis, moreover, 

lim — (x, r) = for any x$ 6 (— oo, oo). 

(x,r)-*(xo,o) r 

More precisely, there exists G (0, 1) such that 

—1 <C(|I 2 -Ji|) 

- r lcP((h,l 2 )x(0,h )) 
holds for < ho < j and any real numbers l\ <l 2 - 



(59) 



(60) 



(61) 



Proof: Since ip £ C 2 '^{Q) f| C 1 '^ U T 2 ) satisfies the equation in ([53]) and admits the 
bound (|59p . therefore, the axially symmetric potential 



tp(x,r) = \H 
J(o,f(o)) \ 



—dx - \ H 
r \ 



-i 



— ar 

r 



(62) 



is well-defined and path independent except on the symmetry axis {r = 0}. Moreover, 



(f x = [ H 



-1 



<Pr 



H 



fx 

r 



(63) 



Therefore, by ([59 



\V<p\ < C in O = {0 < r < f(x)}, 



(64) 



which implies that ip can be extended to as 



ip(xo,0) = lim ip(x, r), Vxo G (— cxd,oo). 

(x,r)£U,(x,r)^(xo,0) 
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Note that the axially symmetric potential (p induces a 3-D potential function 



${x,y,z) = <p(x,Vy 2 + z 2 ) (65) 



which is defined on the three dimensional domain D. Then, 

y z 

$ x = tp x , $ y = <p r -, $ 2 = <p r - in D , (66) 

and that $ G W^(D)f]C 2 ^(D ). Moreover, it follows from §3§, ([66]), and §7$ that for 
any e > 0, $ solves the equation 

div(#(|V$| 2 )V$) = (67) 



in the three dimensional domain D e = {(x,y,z)\ — oo < x < oo,e < a/ 'y 2 + z 2 < f(x)}. 
Therefore, any rj G Cg°(D), one has 

<?(|V<I>| 2 )V$ • V-qdxdydz 



D 



ff(|V$r)V$ -Vrjdxdydz + / // <7(|V$r)V$ • Vndxdydz 

D e JJJd\D e 



II g(\V$\ 2 )^-r]dS + HI <?(|V$| 2 )V$ -Vqdxdydz 
JJdD £ on JJJ D \D e 

/oo p27v q poo pe P27T 

/ 5(|V$| 2 )— <S>(x,£,6)r]sdxd6 + / / / 3(|V$| 2 )V<I> • Vrjrdxdrd9 
-oo Jo or J_ OQ Jo Jo 



as e — ► since V<3? is bounded. Therefore, 

<?(|V$| 2 )V$ • Vrjdxdydz = V// G C^°(D). (68) 



Thus, $ is a weak solution of equation (|67j) in D. Since (|67j ) is elliptic due to (|38|) . thus the 
standard elliptic regularity theory, [ID] , shows 



$ G C°°p). 

Moreover, for /c = 1,2,3, satisfies the equation 

Qi ((5(|V$| 2 )<% + 2~g'(\V^\ 2 )d^d j ^)d J (d k ^)) = 0, (69) 
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which is uniformly elliptic due to (|38|) . Thus, by Nash-Moser iteration, there exists a 
f}\ G (0, 1) such that for suitably small positive constant h, 



Therefore, 



[ d k$} C Pi(B h ((x,0,0))) ^ C\\V$\\L<x>(B 2h ((x,0,0)))' 



\dy$(x, y, z) - d y $(x, 0,0)1 < C(y 2 + z 2 f l/2 for r = {y 2 + z 2 ) 1 ' 2 < h. 



It follows from (1661) that 



\ip r 



(x,r)cos6- $ v (x, 0,0)| < Cr^ 1 , for all 9 £ [0,2vr). 



Thus, 



Similarly, §> z (x,0,0) =0. Therefore 



$ y (x,0,0) =0. 



\tp r (x,r)\ < Cr Pl . 



Thus, 



Furthermore, (fTOj) yields 



lim <Pr( x , r ) = 0. 
(x,r)— >(a:o,0) 



[V<p}ch((h,l 2 )x(0,h)) - C (\ l 2 ~ 



So the desired estimates (f60|) and (fUTj) follow. 
This finishes the proof of the Lemma. 



(70) 



(71) 



(72) 



□ 



Step 9. Removal of cutoff. Combining (j5T|) . ([58]) and ([6T|) yields 



< Cm, in fi. 
If m is sufficiently small, then Cm < m, therefore 

w , - 

< m. 
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Consequently, if) solves the problem (|34p . and moreover, which is uniformly subsonic. 

Step 10. Existence of 3-D subsonic flow. It follows from the proof of Lemma [5] and step 
1-9 that there exists a three dimensional subsonic solution to problem (jlip which satisfies 
CC1 and (HU). 



4 Existence of The Critical Incoming Mass Flux 



In this section, it will be shown that there exists a critical value rh such that the flow is 
always subsonic when the three dimensional mass flux rriQ is less than rh. To achieve this 
goal, we first show that the flow approximates to uniform flows at far fields. 
Let 



H(s) 



H(s) if, s < 4 



ifs>p±!) 2 



(73) 



2 ) )i ° ^ \ 2 

be a smooth decreasing function, where sq G (0, 1). It follows from the proof of Theorem Q] 



that there exists a solution ip to the problem 



V> = o, 

if) = m, 

for any m > 0. Moreover, if) satisfies 



div ( H 



r 



in O, 

on Ti , 
on T 2 , 



(74) 



(75) 



if; < Cr\ 

If the wall of the nozzle tends to be flat at far fields, i.e., / satisfies (|15p . then solutions to 
(|74p approximate to uniform flows at far fields, as is described in the following lemma. 

Lemma 5 Suppose that f satisfies Illty) and \15\). Let if) be a solution to |7^| ) and satisfy 
J75[ ), Then for any e > 0, there exists a constant L > such that 

< e, if x < —L, 



V ib . . . 

(x, r) — (0, 2m) 

r 



and 



Vip 



[x,r 



< e, if x > L. 
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Proof: We start with a special case. Assume that f(x) = a if x > Lq. Set if^k(x, r) = ip{x + 
k,r)xux,r)\x>-k+L +i,o<r<a}- It follows from the standard Holder gradient estimates and 
Lemma[3]that there exists ai S (0, 1) such that for any compact set K C (—00, 00) x [0, a], 
IIV ; fc|lc 1 ' a 2(_ft:) < C(K) for k sufficiently large, where C(K) does not depend on k. Therefore, 
by Arzela-Ascoli lemma, there exists a subsequence ip^, — ► in C 1,Q3 (/C) with 03 < ai- 
Moreover, tpQ solves the following boundary value problem 

£-0o = div((£(|^&| 2 ))- 1 ^) =0, in£ = {(z,r)| -cx)<x<oo,0<r <a}, 

ipo(x,a) = m, if — 00 < x < 00, 

ipo(x,0) = 0, if — 00 < x < 00. 



Furthermore, thanks to ([75]) . t/>o satisfies the estimate 

^0 < O 2 in £ - 
In fact, problem ([76]) and ([77]) has a unique solution 

„2 



mr 



(76) 



(77) 



(78) 



which follows from a simple comparison argument in pQ. 

Since the solution to problem ([76]) and ([77]) is unique, therefore, for if = [—2, 2] x [/j, a], 
then 

- -» for V/i > 0, for fi < a 3 . 

r T 0([-2,2]x[ft,a]) 

By the definition of ipk and ([78]) , this is equivalent to 

— ► as A; — > 00 for V/t > for // < 03. 



^-(0,2m/a 2 ) 
r 



C([fc-2,fc+2]x[/i,a]) 

In the general case that the wall of the nozzle is not flat at far fields, one can set 
tpk(x, r) = ip(x + k, r)x{( x ,r)\x>-k+i,o<r<f(x+k)}- Then it follows from a similar analysis that 
ipki ~~ V'o m C 2 ' a3 (K) for any compact set K C (—00, 00) x (0, a), here K may not touch the 
boundary r = a, but ipo still satisfies the same boundary value problem ([76]) and estimate 
([77]) . Therefore, 

V ^ - (0,2m/a 2 ) 

C^([k-2,k+2]x[h,a-h]) 



as k — > 00 for V/t > for /i < 03. 



(79) 
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However, away from the symmetry axis, ip possesses Holder gradient estimates, conse- 
quently, there exists 04 > such that 

"^1 < a (so) 

. r \c a 4({(x,r)\k-2<x<k+2,r>a-h}) 

Following the same argument in Section 3, one can show that there exists /3 2 £ (0, 1/2) such 
that 

'^1 < a (si) 

. r J C 3 2((fc-2,fc+2)x[0,/i]) 
It now follows from estimates (j79j) . (|80[) and (|81jl that the flow approximates to uniform 

flows at far fields. Indeed, by ([80]) and ([8Tj) . 



(x,n) (x,r 2 ) 

r r 



< C(h% 2 +/ig 4 ) if n,r 2 > a-/i or < ri,r 2 < /i - (82) 



Thus, Ve > 0, it follows from (1821) that there exists h > such that 



(x,n) (x,r 2 ) 

r r 



< 



3' 



V0 < ri, r 2 < /i. 



(83) 



and 



(x,r 2 ) 



< 



3' 



Vri, r 2 > a — h. 



On the other hand, there exists L > such that 



(84) 



— (x,r) - (0,2m/a 2 ) 
r 



eh h 
-3' 2 <^<«-2^> L - 



(85) 



Thus, combining (j83j) . (jSij) and ((BSD y ields 



x,r - 0,— 
r a z 



< e, Vx > L. 



Similarly, we have 



(x, r) — (0, 2m) 



< e, Vx < -L. 



This implies that the flow approximates to uniform flows at far fields. 



□ 



With the help of Lemma [SJ one can show the uniqueness of uniformly subsonic flows. 
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Lemma 6 Suppose that f satisfies lllty) and t!5\) . then uniformly subsonic flows to problem 
(34\ ) are unique. 

Proof: The proof is quite similar to the proof in [13]. Suppose there are two uniformly 
subsonic flows ipi and tp2 which satisfy 

r 

Since ipi (i = 1, 2) satisfies the equation 

*(MW))>)-° 

where H is defined in ([75]) . It is easy to check that tp = ipi — ip2 satisfies an equation of the 
form 

Ctjj = aij(x,r)dijip + bi(x,r)diip = 0. 

By Lemma O the flows corresponding to ipi and ip2 approximate to same uniform flows at 
the far fields, therefore, for any e > 0, there exists a L > such that r)| < e if \x\ > L. 
Thus by maximum principle, |^| < e, Ve > 0, since tp = on T\ and T2. Since e is arbitrary, 
so V? = 0. 

This finishes the proof of the Lemma. □ 

With the help of Lemma Lemma El and going back to the original three dimensional 
flows, we can show in the same way as in |13| that there exists m such that as mo — > rh, 
M(mo) — ► 1. Furthermore, it follows from the comparison principle by Gilbarg[7] that as 
mo 1 m, M(mo) T 1. 

5 Properties of Subsonic Flows 

In this section, as same as the case for plane flows, we will obtain some properties of subsonic 
axially symmetric flows, which are useful to show the existence of subsonic-sonic flows. 



< s < 1. 
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It follows from the proof of Lemma E] that for axially symmetric subsonic flows, problem 
(fTTj) and (fMJ) are equivalent. Thus, in this section, we will use two descriptions simultane- 
ously. 

First of all, as for the plane flows in [13], the axial velocity is always positive. 



Lemma 7 Suppose that f satisfies and M5\) . Let ip be a uniformly subsonic solution 
to problem {3$ , then 

u > in D. 

Proof: Since the flow is uniformly subsonic, one can assume that 



sup 

n 



< s < 1. 



Define g as in (|37p with the help of H in (|73p . and <!> as in (|65p . Then $ satisfies 
div(5(|V$| 2 )V$) = a lj (V^)d ij ^> = 0, in D. 

Set u = & x . Then 

dijCVtydijU + Dp^ijiV^dij^dkU = 0, 

which is a uniform elliptic equation of u. Since ip = m on the solid boundary r = f(x), 
therefore, ip x + ip r f'(x) = 0. On the other hand, ip attains its maximum on the solid 



dtp 

domain Q. Since 



boundary, therefore, by Hopf lemma, J| > 0, where N is unit outward normal to the 2-D 



thus, ^ r > 0. By the definition of ip and $, *0 r > is equivalent to <3? x > on the 



three dimensional solid boundary {(x, r, z)\f(x) = \Jy 2 + z 2 }. Since the flow approximates 
to uniform flows at far fields, moreover, $ x — > {G~ 1 (^r)} 1 / 2 as x — > oo, and $ x — > 
{G _1 (4m 2 )} 1//2 as x — > — oo,, therefore, by the maximum principle 



u > in D. 



22 



Therefore, the proof of the Lemma is complete. 



□ 



Since the axial velocity is positive, then we can define the flow angle by 

V 

lo = arctan — . (86) 
Moreover, we have the following estimate on the flow angles. 

Lemma 8 Suppose that f satisfies Illty) and M5\) . Let ip be a uniformly subsonic solution 
to the problem (3$), then the angle lo defined by KEfy) satisfies 



u<u<u>, (87) 

where 

lo = min{0, inf arctan _f'(x)}, u> = max{0, sup arctan f'(x)}. 

X X 

Proof: The basic idea for the proof of the lemma is the same as that in [13], i.e., using 
hodograph transformation to obtain an elliptic equation for the angle, then the estimate 
([871) will be obtained by a comparison principle for elliptic equations. However, for axially 
symmetric flow, this procedure is more involved. 

Let us first go back to the equations for axially symmetric flows (see (plj) and (f25|) ), 
which reads 

{rgU) x + (rgV) r = 0, 

U r - V x = 0, 

where U = y^, V = — j^, g = g(q 2 ), q = VU 2 + V 2 , and g is defined by ([5]). Since tp 
satisfies (fM|) . so ip as in ([62]) is well-defined, moreover, 

o{x, r) 

and which is strictly positive for r > 0. For r > 0, the mapping (x,r) i— » ((p,ip) is a local 
differmorphism. In fact, the mapping is globally invertible. Indeed, suppose that there are 
two points (xi,r%) and (£2,^2) such that ip(xi,r\) = ^(x2,r2) and ijj(xi,ri) = ip(x2, r^)- If 
ip(x\, ri) = ip(x2,r2) = 0, then it is obvious that xi = X2 and r\ = r2 = due to maximum 
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principle and Lemma UJ Let ■0( x i> r i) = V'( 2; 2,^2) = d > 0, then (x\,ri) and (x2,T2) are 
both on the streamline defined as follows 

f s =U{x,r), 
f s =V{x,r), 
x(0) = x\, r(0) = 

Moreover, this streamline is uniformly away from the symmetry axis, then it follows from 
the argument in [13J that x\ = X2 and r\ = T2- 
Now direct calculations show that 



r{(g{q 2 )U) x + (g{q 2 )V) r ) + g{q 2 )V = r qg(q 2 )(l 



+ rg(q 2 ) 2 q 2 uj^j + g(q 2 )q sin uj, 



U r -V x 



-q 2 0J v + rg(q 2 )qq 4 ,, 



where c is the sound speed. Therefore 



rg(q 2 ) 



which can be rewritten as 



rg(q 2 



1-2 



,2 ^Ip 



+ 



r(l-S) 



/ g \ . , rg(q 2 )q , 



+ 



cos a; sin a; 
-wj, — tv, = 0. 

r(l-g) r*(l-g) * 



It follows from definitions of 93 and ^ that 



Substituting flUD into (JBSJ) yields 



ifx-ipr - ip x (f r rg(q 2 )q 2 ' 



1 1 \ , ( r g(q 2 )q 



+ 



cosw 



1 d 

r dq ' 



. ,. 2 . 2 . sma> 
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■ sinu = 0. 



(89) 



(90) 



(91) 



2 

Since the flow is subsonic, 1 — % > 0, therefore, equation (f9Tj) is an elliptic equation. Note 
that, by Lemma cj G (— §,§)• Moreover, in the domain Q + = {to > 0}f]U, it follows 
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from Lemma [7] that uj satisfies that 

n ( Q \ . frg(q 2 )q \ Id 2 , 2n sinw cosu ^ n 

W? 2 ) V<p A rdq rg(q 2 ) v r(l - g) V 

By the maximum principle, Theorem 3.1 in [10J, 

supu < sup u = maxjsupcj, 0}. (92) 
Similarly, in the domain = {uj < 0} f] £1, uj satisfies 



q \ ( rg{q 2 )q \ 1 d g 2 ^sina> cosw 
by the maximum principle, 



/ u \ 1 1 y\u iq \ 1 a . q amuj <ju» uj 

rg{q) J v ), rd 1 c r r(l-|r) 



inf uj > inf uj = minjinf uj, 0}. (93) 

n- wean- 

Combining (|92l) and (f93l) together, we have 

mm(mfa>, 0) < u < max(supu;, 0). 
9^ an 

This finishes the proof of the Lemma. □ 



At the end of this section, we would like to study the relationship between flow speed 
and incoming mass flux. 

Lemma 9 (a) Let < m± < 1712 < rh. Suppose that tpi are uniform subsonic solutions 
to associated with the incoming mass flux rrii(i = 1,2). Then 

|VVi(x,/(:e)| < |VV 2 (x,/(x))|, e R 1 (94) 

(b ) Both the supremum of flow speed and infimum of horizontal velocity for uniformly 
subsonic solutions to {3$ are monotone increasing with respect to the incoming mass 
flux. Moreover, for any given m E (0, rh), there exist a positive constant 5 = 8(m) > 0, 
such that if m £ \m,rh), then 

q(m) = inf I I > <J_ 
n r 
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Proof: 

(a) It is essentially proved in [7J, however, regularity of solutions in our case near solid 
boundary is only C 1,Q , we can not use Hopf Lemma for derivative of solutions. 

Set tp = tp2 — ipi , then ip satisfies 

Aij(x,r)dijip + Bi(x,r)diip = 0, 

where A^j, Bi £ (£1). Therefore ip achieves its maximum either on the solid 
boundary or at far fields. Note that ip tends to uniform flow at far fields, thus ip 
achieves its maximum on the nozzle wall. Therefore, 

on 

where n is the unit outer normal of the nozzle wall. On the other hand, ipi achieve 
their maximum on the whole nozzle wall where they are constants, so 

dipi dipi 

— — > 0, and — — = 0, on T2, tor 1 = 1, 2, 
on or 

where r is the tangential direction of T2. Thus 

■ V^2 2 > V^1. 2 
r ~ r 

(b) It follows from the proof of Lemma El for any m S (0, m), there exists 5 = 5(m) such 
that the corresponding solution ip with m satisfies 

dib 

inf h= >S. 
n or 

Furthermore, ^ achieves its infimum either on the nozzle wall or at far fields. 
Note that on the nozzle boundary 

dijj dij) 



8„ Or^^- 



therefore, for any m G (m, rh). 



> ^- > mf ^- > mf ^= > 5 

r or or or 
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□ 



Collecting all these lemmas together, we finish the proof of Theorem [2] 

6 Subsonic- Sonic Flows 

To take limit for mo — ► m, let us recall the compensated compactness framework in [13J. 



Theorem 10 Let w £ (x,r) = (q £ , co £ )(x, r) be a sequence of functions satisfying the Condi- 
tions (C): 

(C.l) < 5 < q £ (x,r) < 1 a.e. in Q for some positive constant 5. 
(C.2) \uj £ {x,r)\ < uj < |, for some constant uj independent of e. 

( C.3) d x rj±(w £ ) + d r A±(w £ ) are confined in a compact set in H7 c (0) for the momentum 
entropy- entropy flux pair 

(r/ + , A + ) = (pq 2 cos 2 uj + p, pq 2 sin to cos to), (rj-, A_) = (pq 2 sin oj cos uj, pq 2 sin 2 ui + p), 

where p = p{p), and p = g(q 2 ) is determined by ^ through Bernoulli's law. 

Then there exists a subsequence {w £k } of {w 6 } and w(x,r) = (q,u>)(x,r) such that 

(q £ \uj £ «) -> (q,uj), (95) 
q £k cosu) £k — > qcosu), q £k sina; efc — > gsinu;, (96) 
g((q £k ) 2 )q £k cos oj £k — > g(q 2 )qcosoj, g{{q £k ) 2 )(f k sin oj £k —>■ g(q 2 )q sinoi, (97) 



where all the convergence in 195\) -i 91) are almost everywhere convergence, and w = (q,uj) 
satisfies 

< 5 < q(x,r) < 1, 
|w(a;,r)| < uj. 

Remark 6 The strong convergence of velocity field (U,V) = (qcosu, qsinuj) instead of 
(q,uj) was first proved in [2j. Since we have good control on flow speed, we can also get 
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strong convergence on flow angles. The difference between assumptions on Theorem 1 1 1 and 
Theorem 1 in [2] is that they use one more entropy-entropy flux pair instead of the condition 
on the lower bound on flow speed. 



Let i\) satisfies (|3lj) . Set, as before, 

U=^, V=-^, q 2 = U 2 + V 2 , p = g(q 2 ). 
By direct calculations, it is easy to see that (p, U, V) satisfies 

( P U 2 +p(p)) x + (pUV) r = -^., 

(pUV) x + (pV 2 +p(p)) r --^ 



r 



Let 777.Q —* rh, and (q £ ,uj £ ) be the solutions to (|34p corresponding to 7«q, then away from 
the symmetry axis, on any compact subset in $7, (— SiilAJLL- — ; _ g((g Llff ) ) are uniformly 
bounded, therefore, precompact in H^. Furthermore, due to Lemma [SJ the associated flow 
angles satisfy the estimates 

7T 7T 

— — < min(inf arctan / (x),0) < u £ < max(sup arctan / (x),0) < — . 

Therefore, conditions (C.l), (C.2) and (C.3) in Theorem [10] are all satisfied. Hence it follows 
from Theorem 1101 that (U £ , V s ) has weak-* limit (U, V) such that 

(rg(q 2 )U) x + (rg(q 2 )V) r = 0, 

(98) 

u r - V x = 

holds in $7 in the sense of distribution, where q 2 = U 2 + V 2 , and g is defined in ([5|). 

We now verify that (p, U, V) gives rise to a global subsonic-sonic weak solution to (fTTT) 
on D. Thus, define 

P(x,y,z) = g(q 2 )(x,r), u(x,y,z) = U(x,r), v(x, y, z) = V(x, r)-, w(x, y, z) = V{x, r)-, 
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where r = \Jy 2 + z 2 . First, note that for 7] G Cq°(Dq), 



urjy — vrj x dxdydz 

D 

2n poo pf(x) f smi 

U(x, r)(r] r (x, r, 9) cos 6 — r]e(x, r, 6) ) I rdrdxdO 



J-ooJO 

2tt poo pf{%) 



V(x, r) cos Or] x (x, r, 6)rdrdxd0 

J-ooJO 

/ / r(U(x, r)r] r cos 9 — Vrj x cos 6) — Ur]g sin QdOdrdx 
-oo Jo Jo 

00 rf( x ) / Q r 2lT Q r 27T 



oo JO 



( d f 27T d f 2n \ 

\U(x,r)— / rrj(x, r, 6) cos 6d6 — V(x, r)— / rr](x, r, 6) cos 6d6 ) dxdr 
V or J dx J ) 



0, 



where we have used the facts that 

(■2tt 



rr)(x,r,9) cos 0d6 G Cg°(fi) 



o 



and (|98p holds in the sense of distribution in f2. 
Define C = C( s ) such that 

C(a) = 1 if s > 1, C(s) = 0, if « < ^ K'l < 2 > and (gC c 

and set Q(s) = £(§), where < 6 < b. 

Then for any rj G Cg°(Z)), Cs(r)r/ G Cg°(-Do), thus 



itr^ — vr] x dxdydz 
u(Cs(r)rf) y - v(Cs(r)rj) x dxdydz 



- (s(r))r])y - v((l - Cs(r))v)xdxdydz 

D 

2n poo pS 

d9 / (^(l-C«(r)) tf »7 + ^(l-C«(r))»fe-^cose(l-C^))»fe)r<irdr 

J-ooJO 

-» 

as 5 -> 0. Thus for any rj G Cg°(D) 

itryy — vrj x dxdydz = 0. 
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Similarly, one can show that 

j j j UT] Z — wq x dxdydz = 0, 

and 

j j j vq z — wr]ydxdydz = 0, 

for any 77 G Cq°(D). 

Finally, we check the continuity equation. For any rj £ C£°(D), 

pui]x + pvrjy + pwr] z dxdydz 

D 

pu((s(r)r]) x + pv(Q(r)r])y + pw(Cs(r)v)zdxdydz 

D 

+ j ' J r J pud 1 ~ Cs(r))f]) x + pv((l - (s(r))rj) y + pw({l - Cs{r))r]) z dxdydz 
Note that (ED shows that 

P £ u £ ((sv)x + P £ v £ ((5v)y + p £ w £ (( S 7]) z dxdydz = 0, 

J J JD 

while d9TD yields 

lim j j! Cs(f(x) - r) (p £ u £ ((sv)x + p e v £ ((sv)y + P £ w £ ((srj) z ) dxdydz 
= HI Csifix) - r) (pu(C s (r)rj) x + pv(C s (r)rj) y + pw(( s (r)ri) z ) dxdydz, 

where (p £ , u £ ,v £ ,w £ ) denotes the three dimensional flow associated with the incoming mass 
flux m £ . Therefore, 



pu(( s {r)ri) x + pv((s{r)ri)y + pw((s(r)r)) z dxdydz 

D 

1 - Csifix) - r))(pu{(sir)v)x + pviCsir)v)y + pw (C<5 (r)) rj) z dxdydz 

D 



D 

+ lim j II (sifix) - r) ip £ u £ iCsv)x + p £ v £ iCsv)y + P £ w £ (Csv)z) dxdydz 




+ I IJ Csifix) - r)ipu(( s (r)ri) x + pv((s(r)v)v + pwiCsir))v)zdxdydz 

(! - Csifix) ~ r))ipuiCsir)v)x + pviCsir)v)y + pwiCsir))r]) z dxdydz 
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D 



(1 - Cs(f(x) - r))(pu((g(r)rj) x + pv(Cs(r)rj) y + pw{Cs{r))r]) z dxdydz 



+ lim / / 1 (p £ u £ (Csv)x + P £ v £ ((sri)y + p £ w £ {(sr]) z ) dxdydz 



£^0 



+ lim / / / (ft (/(a;) ~r)-l) (p £ u £ {(sv)x + p £ v £ (Csv)y + p £ w £ ((sv)z) dxdydz 



y 

D 




(1 - Cs(f(x) - r)){pu((s{r)rj) x + pv(( 5 (r)r)) y + pw(( s (r))r]) z dxdydz 

D 

+ lhnj I J (&(/(*) - r) - 1) (p £ n £ (C 5 r/)x. + p £ v £ (C^)s/ + /^ £ (C5r?) 2 ) dxdydar. 



Thus 



ptiTya; + pvrjy + pwq z dxdydz 

- C^)) 7 ?)^ + /Of ((1 - Cs{r))r]) y + - ( s (r))r]) z dxdydz 

+ / / / (1 - C<5(/0*0 - r )) (pu((s(r)v)x + pv(( s (r)ri)y + pw(Q(r)r]) z ) dxdydz 



D 



D 



D 

+ lim y jj (Uifix) - r) ~ 1) (jftfiCM)* + P £ v £ ((sv)y + P £ w £ ((sv)z) 




as 5 — > 0. Therefore, for any test function 77 E C%°(D), 

purjx + /w% + pwrj z dxdydz = 0. (99) 



Moreover, the equation (|99p also implies that (u,v,w) satisfies the boundary condition ([9]) 
actually as the normal trace of divergence measure field (pu, pv, pw) on the boundary in the 
sense of Anzellotti[2j. So we finish the proof of Theorem 3. 

Further characterizations of the subsonic-sonic flow we obtained are left for future. 
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